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Structure and production of doubly charmed tetraquarks T cc (ccud) are studied from the 
viewpoint of color configurations. Based on the diquark correlation, the tetraquark T cc with 
I(J P ) = 0(1 + ) is considered to be stable against strong decay. We discuss that the mixing prob- 
ability of color antitriplet and sextet cc components in T cc is suppressed by l/m^, so the two 
configurations are separately realized in the heavy quark limit. Utilizing the nonrelativistic QCD 
framework, we evaluate the production cross sections of T cc in electron-positron collisions. The 
momentum dependence of the cross section of color antitriplet is found to be different from that 
of sextet, which can be used to discriminate the color structure of the T cc states in experimental 
measurements. 

PACS numbers: 14.40.Rt, 14.40. Lb, 13. 66. Be, 12.39.Jh 



Study of exotic hadrons is one of the most interesting 
subjects in the quark and hadron physics in these years. 
The experimental discoveries of the exotic hadrons, such 
as X, Y, Z ± in the charm sector and Yb, Zj^ in the bot- 
tom sector, have motivated many researchers towards the 
study of those heavy- flavor exotic hadrons [THl] . The ex- 
otic hadrons have many unexpected properties, such as 
masses, decay widths, branching rations, which are hard 
to be explained in the conventional quark model. The 
doubly charmed (C = 2) tetraquark T cc with four quark 
configuration ccud (two charm quarks c's and light up 
and down antiquarks u and d) is a new interesting can- 
didate, because it is a genuine tetraquark hadronF] The 
mass of T cc as a compact tetraquark state has been the- 
oretically studied from the quark model [6- 29J and from 
the QCD sum rules [3U]. It is studied also as molecule-like 
DD, DD* and D*D* bound or resonant states in hadronic 
molecule picture [3TH35] , One of the interesting proper- 
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1 in literature, another definition of T cc as ccud has also been 
used. Although we deal with the ccud state in this paper, all the 
discussions can be equally applied to ccud state. 



ties of T cc is its color configurations. A quark-quark pair 
in T cc can be combined, not only in color antitriplet (3) 
configuration, but also in color sextet (6) configuration, 
the latter of which does not exist in normal hadrons with 
qqq or qq. If T cc exists, it helps us to understand the 
quark-quark interaction in the color channel which is not 
accessible in normal hadrons. To research the existence 
of T cc , it is required to search for T cc by experiments in 
high energy accelerator facilities [33 [37] . In the present 
study, we consider T cc as a compact tetraquark state and 
study how T cc can be produced in electron-positron col- 
lisions. 

In literature, there have been many quark model calcu- 
lations of the mass of T cc . In shortly summarized, their 
results can be understood by the diquark picture [55] 
where diquarks are compact objects. The attraction for 
the diquark is supplied from the color-spin interaction, 

ff «- = E^j(-|)Wi^- w 

Here Ch is the coupling strength, m t is the mass of the 
quark i, A; is the Gell-Mann matrices operating to the 
color of the quark i, and Sj = <x;/2 is the spin operator 
with the Pauli matrices <?, operating to the spin of the 
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quark i. The factor —3/8 is multiplied for normaliza- 
tion (— 3/8\i-\j = 1) for the pair of quarks i and j in 
color antitriplet (3) channel. The values of Xi-Xj Si-Sj are 
dependent on the color-spin channels as summarized in 
Table [j] We note C H is given as Ch = v (S(rij)), where 
vq is related to the interaction constant, rij is the dis- 
tance between the quarks i and j, and the expectation 
value of the delta function is given from the wave func- 
tion of the relative coordinate for the quarks i and j. 
The parameters are given by Ch = Cb (2Ch = Cm) for 
quark-quark (quark-antiquark) pair with Ch/to^ = 193 
MeV (CmK = 635 MeV ), m u = m d = 300 MeV 
and m c = 1500 MeV to reproduce the masses of known 
hadrons [201 [35] • Note that Cm is larger than Cb, re- 
flecting the color factor by two and the difference of the 
quark wave function. For example, in A c , there is an at- 
tractive ud diquark (''good" diquark) with color 3, spin 
1 So in isospin I = 0|j In E c and £*, there is a repulsive 
ud diquark ("bad" diquark) with color 3, but spin 3 Si in 
isospin 1=1. 

As a remark, we note that Cb for heavy-heavy (cc) 
quark pair should be different from those of light-light 
and light-heavy quark pairs. This is because the wave 
function of the heavy-heavy quark pair is compact and 
its value at the center of mass is increased. Hence the 
expectation value (S(r)) by the wave function (r the dis- 
tance between the two quarks) is increased. Therefore, 
Ch for charm quark pairs should be different from the 
other cases as Ch = C cc with C cc /m^ = 39 MeV, which 
causes only the minor difference [20J [25] . 

Now let us see how T cc can be a stable particle in 
the diquark model [201 l2"5] . The flavor dependence of 
the color-spin interaction is important to the stability of 
T cc . We note that the color-spin interaction H{ nt is pro- 
portional to the quark mass factor 1/mjWj. It indicates 
that the light quarks are affected much by the color-spin 
interaction, while the heavy quarks are not. We thus 
consider a state with "good" ud diquark and cc diquark 
with color 3 and spin 3 Si. The repulsion of cc diquark 
is very small due to the suppression with the 1/m;: fac- 
tor, while the ud diquark induces a strong attraction to 



TABLE I. The expectation values of (— 3/8)Ai • Aj ■ Si ■ Sj for 
quarks i and j with spin s = 0, 1 and color 3, 6. 
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2 It is important to note that the attractive ud diquark is a source 
to form the color superconductivity in quark matter at high den- 
sity l38H40j . 



make a bound state. In this picture, the quantum num- 
ber of the ground state of T cc is uniquely determined as 
I(J P ) = 0(1 + ). The lowest two-meson threshold in this 
channel is the D and D* mesons in s wave. We can esti- 
mate the binding energy of T cc from the DD* threshold 
as 

B.E. = (4— + \—) (A— 2 + 7%) 

\ 4 m u m c Am u m c J \ 4 m; 4 ra~ J 

~ 71 MeV, 

where the first term is the sum of the energies from color- 
spin interaction in D and D* mesons and the second term 
is the energy from color-spin interaction in T cc . The 
above result means that the mass of T cc is 71 MeV below 
the DD* threshold. Thus, T cc can be a stable bound state 
as a genuine flavor exotic hadron. 

We consider a new state of T cc whose color configu- 
ration is different from the conventional T cc discussed 
above. The new state of T cc has the configuration of cc 
(ud) pair in color 6 (6) and spin 1 So ( 3 Si). We find that 
the ud pair in color 6 and spin 3 Si in isospin I = is 
attractive, although its strength is smaller (factor —1/8 
from Table |l| than that of the ud diquark in color 3, 
spin 1 So and isospin 1 = (factor —3/4). Denoting the 
quantum numbers of the cc pair, we shall refer to the 
conventional (new) state as T cc [3, 3 Si] (T cc [6, %]). The 
quantum number of T cc [6, 1 So] is again I{J P ) = 0(1 + ). 

In principle, T cc [3, 3 Si] and T cc [6, -"-So] can be mixed 
because both of them have the same quantum numbers 
I{J P ) = 0(1+). However, the spins of cc in T cc [3, 3 Si] 
and T cc [6, 1 So] are different, so the transition requires 
the spin-flip of a heavy quark which is suppressed by 
l/m c . Therefore, the mixing probability of T cc [3, 3 Si] 
and T cc [6, 1 So] is suppressed by order of 1/m;! and we ne- 
glect the mixing effect in this work. We can estimate the 
mass splitting between T cc [3, 3 Si] and T cc [6, 1 So] from 
the color-spin interaction iJ; nt as 

M(T CC [6, ^o]) - Af(T cc [3, 3 Sr]) 

= (A^ + _ Cb 1 Ccc\ 

\ 8 m% 8 m% J \ 4 Ami J 
~ 125 MeV. 

Therefore, T cc [6, %] is an excited state of T cc [3, 3 Si]. 
We should remark on the stability of T^G/So]. The 
decay of T cc [6, 1 So] into T cc [3, 3 Si] requires the two-pion 
emission at least from the isospin conservation, but this 
is kinematically forbidden due to the mass splitting be- 
tween T cc [3, 3 Si] and T cc [6, 1 So] as estimated abovej 3 ] 
The mass splitting indicates that the mass of T cc [6, 1 So] 
can lie at 54 MeV above DD* threshold, when the mass 
of T C J3, 3 Si] lies at 71 MeV below DD* threshold (see 
Fig. fl]). Therefore it should be kept in mind that 



The photon emission in electromagnetic interaction gives a small 
width. 
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FIG. 1. Masses of T cc [3, 3 Si] and T cc [6, 1 S ] (solid lines) from 
the diquark model. The DD* and D*D* thresholds are indi- 
cated by the dotted lines. The color and spin configurations 
of quarks are also shown. 



T cc [6, 1 So] may decay to D and D* mesons via strong 
interaction in s-wave fall-apart process. However, the 
decay to DD* requires the color recombination from the 
color sextet cc and antisextet iid diquarks to color sin- 
glet cu and cd mesons, and the size of cc in T cc [6, 1 So] 
is compact because of the heavy reduced mass. These 
qualitative arguments suggest that T cc [6, -"-So] can be a 
narrow state. 

Now let us consider the production of T cc in inclusive 
processes e + e~ — > T cc +X. We use the framework of non- 
relativistic QCD (NRQCD) @TJ [42] which is an effective 
field theory based on the expansion in terms of the heavy 
quark velocity v. The production cross section is factor- 
ized into the short-distance perturbative amplitude and 
the nonperturbative matrix element of the NRQCD oper- 
ators. The framework has been applied to double-charm 
productions with charmonia in the final states |4"3Tl4l)] . 

Here we assume that the factorization is also valid for 
T cc , following the strategy of Refs. [47l [48] where the 
production of the doubly-charmed baryon S cc is studied 
in NRQCD. The cross section is then decomposed as 

da a (e+c- ^T cc [a]+X) 

= ^c1»t(c+c- -> [cc} k a + c + c)(O k (T cc [a])), 
ft 

where da represents the short-distance part of the process 
e + e~ — > [cc]*; + 5 + 5 with the cc pair being projected 
onto the definite color-spin state labeled by a — [3, 3 Si] 
or [6,%] for T cc [3, 3 Si] or T cc [6, %], respectively, and 
k specifies NRQCD operators which are sorted out by 
velocity expansion. Feynman diagrams for the leading 
order inperturbative QCD contributions to da are shown 
in Fig. V2\ The matrix element (O fc (T cc [a])) represents 
the long-distance nonperturbative process of the cc pair 
into T cc . In the leading order of NRQCD, the matrix 
element for each channel is given by a single constant as 



(O k (T cc [a})) 



= U[3,3 Sl ] fora=[3, 3 Si] 
fc=Lo 1 /i[a,is ] for a = I 6 ' ls o] 
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FIG. 2. Feynman diagrams for the leading order contribu- 
tions to the perturbative d<r(e + e - — > [cc]* + c + c) process. 



As we have discussed, since the different color-spin states 
do not mix with each other in the heavy quark limit, we 
calculate their productions separately^ In fact, the color 
6 (3) cc component in T cc [3, 1 S ] (T cc [6, %] ) is included 
in higher order in NRQCD expansion. 

To evaluate da, we assign kinematic variables as 
e + (pi)e~(p 2 ) -> [cc]*(p) + c(p 3 ) + c(p 4 ). We work in 
the center-of-mass frame of the e + e~ collisions with the 
z axis in the beam direction. Defining the x axis so that 
the three-momentum of the produced T cc lies in the xz 
plane, we write the four-momentum of T cc as 

p M =(E p , psin0,O,pcos6), 

where E p = 4m^ + p 2 , p = \p\ is the magnitude of the 
three-momentum of T cc and is the angle of p from the 
z axis. The momenta of c are specified by P3 = — p/2 + q 
and p 4 = — p/2 — q, and q is expressed in the cylindrical 
polar coordinates as q = (qsin6,q y ,qcos0). 

The final expression of the differential cross section is 



da n 



1 



P 



d(9 



dpdcosO (27r) 4 16m c sE p j 

yj A - Bq 2 + Cq 2 cos 2 6' 

A 

B-Ccos 2 B r 



dq 



_ q\M a \ 2 h a 



q y (E 3 + E 4 )' 



q v =- 



4 The final state interaction of T cc [6, ^o] with other hadrons may 
affect on the production of T cc [3, 3 Si] and vice versa. However, 
such a quantitative analysis is not covered in the present study 
and is left for future works. 
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where s is the total energy squared, £^4 = 
V m l + IP3, 4 | 2 , 0' = e-Q, A = js{J- s -2E p ){J- S ~E p )\ 
B = 4(Js - Ep) 2 , and C = Ap 2 . 

The amplitude for the e + e~ — > [cc]^, + c + c process 
A4 a is calculated by the diagrams in Fig. [2j with the 
color-spin projection 

p £l = £ **(§)*%(§) = ^ (2 +-c) ^ (A) ^,, 

P 6 ,m = £ «*(§)*%(§) = ^ +™c) 75C$L„ 

where the sum is taken for each color-spin state, and 
is the polarization vector of T cc , C is the charge conju- 

A S AS 

gation matrix, and $ m 'j.„- = T^ m 'jfc are the normalized 
tensor in color space. 

We fix the total energy at y/s = 10.6 GcV as in Belle 
experiment [49] . Following the discussion of the produc- 
tion of doubly charmed baryon S cc in Ref . [IS] , the other 
constants are chosen to be m c = 1.8 GeV and a s — 0.212. 
In the leading order of NRQCD, the mass of T cc is given 
by Mt cc = 2m c . Note that the mass of charm quark is 
different from that used in the diquark model. We calcu- 
late the differential cross sections normalized by the total 
cross section (l/a a )da a /(dp dcosO) which is indepen- 
dent of the value of the matrix element h a . The results 
are plotted as functions of the T cc momentum p = \p\ in 
Fig. [3]with = 0, 7r/2, and n. Because the amplitude is 
symmetric under the exchange p\ f)p2, the result with 
and that with tt — are identical. In Fig. [3] we see that 
the momentum distributions of the two color states are 
quite different. The T cc [3, 1 So] channel has maximum at 
p ~ 3.5 GeV with short tail in the low momentum re- 
gion, while the distribution of the TccfB, 1 So] state peaks 
at around p ~ 2.7 GeV with appreciable strength in the 
low momentum region. This qualitative difference can 
be used to distinguish the color structure of T cc ; for in- 
stance, the ratio of the cross sections at p = 1 GeV and 
p = 3.5 GeV is quite different in two cases. 

The nonperturbative constants /iraagji and tye^So] 
should in principle be estimated in QCD. In the present 
study, as a simple approach, we use the non-relativistic 

T [3 1 9 

quark model. We assume /irjjssji = \Rcc ' (0)| /47r 

and h[Qis ] = \Rcc° ' (0)| 2 /47r, respectively, where 

Rcc cl ^ Sl] (r) [i?J c cc[6c ' 1Sol (r)] is the radial wave function 
of cc quark pair in tetraquark T cc [3, 3 Si] (T cc [6, 1 So]). To 

obtain the wave function i£H S ' 3Sll (r) and i?J c c c[6 ' 1So] (r), 
we consider the hamiltonian with the harmonic oscillator 
potential J2i<j (~ Tg) ^»"Aj ^j\ 2 where fj is the posi- 

tion of the quark i, and k is the strength parameter of the 
harmonic oscillator potential for quark confinement. The 
color factor Xi-Xj is important to obtain the wave function 
of the cc quark pair, because the wave function should be 
different for each color channel, 3 and 6. Although the 
harmonic oscillator potential is a simple potential in the 
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FIG. 3. (Color online) Differential cross section as functions 
of the magnitude of the three momentum of T cc . Thin (thick) 
lines represent the production of T cc [3, 3 Si] (T cc [6, 1 So]). 



quark model, it will be enough for our purpose when the 
cc wave function contains only the s wave. The masses 
of quarks are set to be the same value used in Eq. |l]). 
The strength parameter k — 0.33 GeV 3 is fixed to repro- 
duce the value of the wave function at the center-of-mass 
for cc quark pair in charmonia, |i? cE (0)| 2 , which is esti- 
mated in the more sophisticated quark model with the 
Cornel-type (Coulomb + linear confinement) potential 
and the spin-spin interaction in Ref. [50] . Here we ob- 
tain |i? cE (0)| 2 = (|^(0)| 2 + 3|i?^(0)| 2 )/4 = (1.18) 2 
GeV 3 with the wave functions |i?^(0)| 2 = (1.39) 2 GeV 3 

and liicg^O)! 2 = (1.10) 2 GeV 3 at the center of mass 
of cc in rj c and 3 ftp, respectively. With all the param- 
eters being fixed, we calculate the cc wave functions in 
T cc [3, 3 Si] and T cc [6, So]- The harmonic oscillator po- 
tential has the property that the cc wave function is ex- 
actly decoupled from the other light quarks (u and d). 
Here we may ignore the color-spin interaction for cc pair, 
because in first the cc wave function is decoupled already 
from the other light quarks, and in second the color-spin 
interaction for the cc quarks is small with suppression 
factor by 1/m 2 . From the harmonic oscillator potential, 
as numerical result, we obtain fya^g^ = 0.089 GeV 3 and 
^[6 c . 1 s ] = 0.054 GeV 3 . By substituting these values, we 
finally obtain the total cross section at Belle energy as 
CT[3 3 Sl ] = 18.8 fb and f[ 6) ig ] = 4.1 fb. We may use 
the other parameter set, for example m c = 1.5 GeV and 
a s = 0.26 from Ref. [H]. Then we find a few times 
larger cross sections cr^ ag^ = 65 fb, <T[ 6) ig ] = 21 fb with 
leaving the momentum dependence in Fig. [3] unchanged 
qualitatively. 

Some remarks are ready. The cross section of the dou- 
bly charmed baryon S cc has been discussed in a simi- 
lar method [57J |3H] ■ When the harmonic oscillator po- 
tential is used, interestingly, the wave function of cc in 
S cc can be estimated as the same with that of cc in the 
tetraquark T cc [3, 3 Si]. Then, we obtain the cross section 
of S cc which is same as that of T cc [3, 3 Si]. The obtained 
S cc total cross section is comparable to the previous re- 
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suits of - 10 2 fb 071 HB]. 

However, there are two issues. First, we have used a 
very simplified quark model having the harmonic oscil- 
lator potential. Although we expect that it will be valid 
to consider only the s-wave state, we will need to ana- 
lyze the wave functions in more realistic potential such 
as the Cornel-type potential for more quantitative dis- 
cussion. Second, we have assumed that the long distance 
quantities hrg sg^ and h[Qis ] as well as hs co for S cc are 
related to the cc wave function. It should be noted that 
the difference of the numbers of the light quarks (two for 
T cc [3, 3 Si] and T cc [6, 1 So], one for S cc ) is not accounted 
for in the present estimate of the cross section. How- 
ever, as the number of light quarks is increased, it should 
be expected that the production of the hadrons becomes 
more suppressed. Thus, the cross section of T cc [3, 3 Si] 
and T cc [6, 1 Sq] should be smaller than that of S cc . For 
more quantitative study, we need to include the light 
quark degrees of freedom in the fragmentation process. 

In summary, we discuss the exotic color configurations 
in doubly charmed tetraquark T cc . We discuss that T cc 
with color 6 configuration does not mix with conventional 
T cc with color 3 in the heavy quark limit. We evalu- 



ate the production cross section of T cc in the electron- 
positron collisions by NRQCD. As a result we find that 
T cc with color 6 has different momentum and angular de- 
pendence from that of T cc with color 3. This study will 
be useful to pin down the color configurations in exotic 
hadrons from experimental measurement. 
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